New Physics and the Proton Radius Problem 
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Background: The recent disagreement between the proton charge radius extracted from Lamb 
shift measurements of muonic and electronic hydrogen invites speculation that new physics may 
be to blame. Several proposals have been made for new particles that account for both the Lamb 
shift and the muon anomalous moment discrepancies. Purpose: We explore the possibility that 
new particles' couplings to the muon can be fine-tuned to account for all experimental constraints. 
Method: We consider two fine-tuned models, the first involving new particles with scalar and 
pseudoscalar couplings, and the second involving new particles with vector and axial couplings. 
The couplings are constrained by the Lamb shift and muon magnetic moments measurements while 
mass constraints are obtained by kaon decay rate data. Results: For the scalar-pseudoscalar model, 
masses between 100 to 200 MeV are not allowed. For the vector model, masses below about 200 MeV 
are not allowed. The strength of the couplings for both models approach that of electrodynamics 
for particle masses of about 2 GeV. Conclusions: New physics with fine tuned couplings may be 
entertained as a possible explanation for the Lamb shift discrepancy. 



I. INTRODUCTION 

The recent measurement of the muonic hydrogen Lamb 
shift [l[ yielded a proton charge radius 5a smaller than 
the 2006 CODATA value available at the time of its pub- 
lication @, and la smaller than the 2010 CODATA up- 
date Q , which incorporates the latest proton radius de- 
terminations from electron scattering [J] . The data used 
in the CODATA determinations is all electronic. A pos- 
sible explanation for the surprising muonic result is that 
overlooked Standard Model processes or new physics are 
at present being wrongly attributed to proton size ef- 
fects. Assuming the previous electronic measurments of 
the proton charge radius are correct, alternative expla- 
nations must lower the muonic Lamb shift by 310 /ieV to 
match the experimental result. 

Several new physics proposals have been considered to 
explain the discrepancy. Jaeckel and Roy as part of 
a larger investigation into deviations of Coulomb's law, 
determined that hidden photons that coupled equally to 
electrons and muons could not explain the discrepancy 
as they would actually cause the proton radius to ap- 
pear smaller in ordinary hydrogen. Tucker-Smith and 
Yavin [(| developed two simple models in which either a 
new scalar or vector particle couples to muons and pro- 
tons. The strength of the particle couplings were set 
by constraints placed by the small discrepancy of the 
muon's anomalous magnetic moment. They showed that 
in first order of nonrelativistic perturbation theory, ex- 
change particle masses of order MeV could produce the 
observed Lamb shift, albeit exchanges with masses this 
light run afoul of neutron scattering data if the neutron 
and proton have similar coupling. Barger et al. Q also 
considered new scalar and vector particles, but suggested 
it would be difficult for them to satisfy additional con- 
straints placed by T, J/tp, n, and ?7-decays. 

Batell et al. § revived the possibility that a hid- 
den photon could be responsible for the Lamb shift dis- 
crepancy by requiring it to couple only to right-handed 



muons. This boson also mixes with the photon so the 
couplings contained additional model dependence not 
seen in other proposals. In order to account for the 
muon anomalous moment constraint, they were forced 
to introduce and fine-tune the mass of a new scalar par- 
ticle. In a second paper Barger et al. @ noted that Batell 
et aVs model does not respect the constraint placed by 
if-decay [l(| if the decay of their hidden photon were 
invisible. 

In this work we explore the possibility that fine-tuned 
particle couplings, free from the phenomenological de- 
mands of hidden photons, can satisfy muon anomalous 
moment and X-decay, as well as other, constraints. We 
consider two separate possibilities. The first contains 
two new particles that interact with muons and protons 
through fine-tuned scalar and pseudoscalar couplings, re- 
spectively. The second contains two new particles that 
interact with muons and protons through fine-tuned po- 
lar and axial vector couplings, respectively. 

Our evaluation of the particle mass and coupling pa- 
rameters proceeds as follows. We begin in Sec. |TT] by 
finding what coupling parameters are needed to obtain 
an extra 310 /xeV muonic hydrogen Lamb shift from the 
exchange of an electrophobic spin-0 or spin-1 particle of 
a given mass. In the nonrelativistic limit, pseudoscalar 
and axial vector particles do not contribute significantly 
to this shift and their couplings remain free parameters. 
Then in Sec. lIIII we confront our models to the constraint 
for the muon anomalous magnetic moment. Polar vec- 
tor and axial vector exchange give opposite contributions 
to the magnetic moment, as do scalar and pseudoscalar 
exchanges, so cancellations can be arranged if the pseu- 
doscalar and axial couplings are tuned. We further con- 
sider in Sec. |IV] the consequences for the decay of K's to 
two or more unobserved particles. This decay is possible 
as a radiative correction to K — > fiv if there is a coupling 
of a light new particle to the muon, and strong experi- 
mental limits are known. In Sec. [V] we make some final 
remarks. 
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II. LAMB SHIFT 

For where there are scalar and pseudoscalar par- 

ticles coupled to a muon and proton, 



(1) 



a coupling much stronger that the Weak coupling, even 
without the (1/4 — sin 2 W ) factors. 

The exchange of a Higgs boson also cannot account 
for the Lamb shift discrepancy. The energy shift due to 
a Higgs boson of mass 125 GeV is roughly eight orders 
of magnitude smaller than what is observed. 



where cf> is a scalar and tp is a pseudoscalar field. The 
potential in the nonrelativistic limit is 
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where M is the mass of the exchanged particle, cj). 

The pseudoscalar contributions are much smaller at 
low momentum transfer. For the 2P-2 S sp litting in hy- 
drogen one gets an energy difference @, [ll| , 
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where (for the muonic case) m r = m^nip/ (m^ + m p ). 

The scalar coupling required to give an extra 310 fxeV 
to the muonic hydrogen 2S-2P Lamb shift, in the case 
Cg = Cg = Cs, is shown by the solid line in Fig. Q] as 
a function of the exchanged mass, up to an exchanged 
mass of 1000 McV. 

For the case of polar vector and axial vector particles 
coupling to the muon and proton, 

C v = - CO^^W - C'a^^'Vi/TsV'm 

- &v4> v i) P i v ii) p - c p A Lp l 'Vyy„7 5 ^p. (4) 
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FIG. 1: The scalar and pseudoscalar couplings needed to sat- 
isfy the experimental constraints. The scalar coupling (solid 
line) is required to give an extra 310 fieV to the muonic hy- 
drogen 2S-2P Lamb shift. The dashed line is the pseudoscalar 
coupling needed to satisfy the constraint placed by the muon 
anomalous moment. We assume the two particle masses are 
identical. 



where 4> v is a polar voector and p u is an axial vector field. 

In the nonrelativistic limit the potential is similar to 
Eq. ©, with vector couplings taking the place of the 
scalar ones and opposite overall sign. The solid line in 
Fig. [5] displays the vector coupling strength as a function 
of the exchanged mass. 

As a side note, these results agree with the known 
Weak interaction or Z-boson exchange contribution to 
the Lamb shift of a given nl state of hydrogen [l2| , 



a{Za) 3 m r 8G F m 2 r fl 
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with the appropriate substitutions. (Just the vector part 
of the Z interaction contributes; further M — > Mz, 
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9v 
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(6) 

and Gf/V2 = g 2 /(8M^ cos 2 Qw)- The numerical result 
for the n = 2 levels is about a million times smaller than 
the extra muonic Lamb shift, speaking to the fact that to 
get an effect relevant to the Lamb shift problem with an 
exchange particle mass in the 90 GeV range would require 
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FIG. 2: The polar and axial vector couplings needed to sat- 
isfy the experimental constraints. The vector coupling (solid 
line) is required to give an extra 310 fieV to the muonic hy- 
drogen 2S-2P Lamb shift. The dashed line is the axial cou- 
pling needed to satisfy the constraint placed by the muon 
anomalous moment. We assume the two particle masses are 
identical. 
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III. MUON MAGNETIC MOMENT 

The muon anomalous moment is accurately measured. 
The theory for the anomalous moment is also quite accu- 
rate, with the bulk of the error coming from uncertainties 
in hadronic contributions. There is a small but persistent 
discrepancy between experiment and theory. In terms of 
0|t = (g - 2) M /2, 

a M (data) = (116 592 089 ± 63) x 1CT 11 [0.5 ppm], 

a M (thy.) = (116 591 840 ± 59) x 10" 11 [0.5 ppm], 

5a^ = (249 ± 87) x 10 -11 [2.1 ppm ± 0.7 ppm]. 

(7) 

The data is from [l3|, E3] and the latest theory number 
is from [loT ]. 

This discrepancy is four orders of magnitude in frac- 
tional terms smaller than the one due to the Lamb shift. 
Every particle that contributes to the Lamb shift also 
contributes to the magnetic moment at the one loop level, 
as in Fig. [3] The contributions of the pseudoscalar and 
axial vector, whose couplings are not constrained by the 
Lamb shift, have opposite sign to those from the scalar 
and polar vector, and can be tuned to respect this much 
smaller discrepancy. 
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FIG. 3: One-loop magnetic moment correction 



For scalar and pseudoscalar particles, we consider their 
masses to be the same. The magnetic moment result is 
known in the literature [ill [TtJ , 
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(The expressions continue nicely to r < 4.) Low and high 
mass limits are 
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Eq. © can be rearranged to solve for (Cp) 2 . The result 
is plotted as the dashed line in Fig. [T] One notices that 
fine tuning must be done to several significant figures at 
higher masses. 

For polar and axial couplings, we also only consider the 
case where their masses are equal. Their contribution to 
the muon's magnetic moment is 
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Rearranging Eq. (fT3)) allows for the evaluation of 
(C^) 2 . The result is plotted as the dashed line in Fig. [2] 
More general combinations of S, P, V, and A are 
also possible. Of note is the model of Ref. Q, which 
involves a vector particle with extra parity violating cou- 
pling to the muon with (in our notation) fixed C v and 
(10) C/t' th&t achieves nne tuning using a scalar, also with 
definite muonic coupling, but with a tunable mass. 
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IV. K DECAY WITH UNOBSERVED 
NEUTRALS 

If a light neutral particle couples to muons, the decay 
K — > fii/cj) is possible, Fig. SJ There has been an exper- 
imental search for multibody decays K — > [iX where X 
contains only neutral particles that are not photons fiol ] . 
The experiment searched for muons in the kinetic en- 
ergy range 60 to 100 MeV — for reference, the muon from 
K -> fw has = 152 MeV or E, = 258 MeV— and 
found a strong limit. To give the result with some pre- 
cision, the experimenters give their detector efficiency 
function D(E,), which is the relative efficiency to de- 
tect an energy E, muon compared to a 258 MeV muon, 
and which is zero outside the stated limits and smoothly 
varying in between. 



K(k) 




FIG. 4: Kaon decay with an extra neutral scalar, <f>, either 
scalar or pseudoscalar. 



For a given decay spectrum dT(K — > /j,X)/dE fl , the 
experimental limit is quoted as [lpl ] 



1 



dT(K -> fjJC) 
dE, 



D{E tl ) dE, < 2 x 10" 



The simple K,2 decay rate is 
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where the kaon decay constant is defined from 

(0| U7 M (1 - 75 )s |0) = V2f K k„ . (20) 
Generically, a three body decay is given by 

T(K fMU(f>) = J dE, dE v l-^l 2 , (21) 

spins 



with integration limits 
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The matrix element for the decay into a muon, neu- 
trino, and particle with both scalar and pseudoscalar 
couplings is 

G F f K V us _ 

x [(Cs - *Cp)Q 2 + ™»{Cs + * C p) H] (1 - 7sM<z) , 

(24) 

where Q 2 = (k — q) 2 = m 2 K — 2m K E l/ . 

The matrix element squared and summed is 
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We evaluate the left-hand-side of Eq. (fT8]) using our con- 
strained Cg and Cp couplings for a given scalar or pseu- 
doscalar mass m^. We note that neither the anomalous 
moment nor the square of the matrix element contain 
terms with both Cg and Cp. Thus, a model with two 
equal mass scalar and pseudoscalar particles is indistin- 
guishable from one which has only a single particle with 
both scalar and pseudoscalar couplings. 

Comparison of our calculated expcrimcntal-cfficicncy- 
weighted decay rate to the experimental limit is shown 
in Fig. [5j A range of scalar masses from about 100 to 
200 MeV is not allowed. For other masses, the couplings 
are not excluded. 

The matrix element for the decay into a muon, neu- 
trino, and a particle with both polar and axial vector 
couplings is, using Cp L = C v ± 

x[C^Q 2 + m^f(}(l- l5 )v(q), (26) 

where e„ is the polarization vector of the new particle. 
This leads to 
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(27) 
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FIG. 5: Mass limits on scalar and pseudoscalar particles due 
to constraints placed by K — > \iX searches. The solid curve 
is the full result, accounting for the experimental efficiency, 
obtained through satisfying the Lamb shift and magnetic mo- 
ment criteria. The contributions of the scalar (dashed curve) 
and pseudoscalar (dash-dotted curve) couplings are indicated 
separately. The experimental limit is the horizontal line, and 
the shaded region is allowed. 



We note that this matrix element squared contains 
terms with CyC 1 ^. Hence one must distinguish the case 
of two particles of equal mass from one parity violating 
particle with polar and axial couplings. The anomalous 
moment has no cross terms, so both the one and two par- 
ticle cases contain the same couplings. However, in the 
decay process, the parity violating case has the possibil- 
ity of destructive interference. 

The results for both cases are shown in Fig. [6l us- 
ing couplings obtained from the fine tuning of the muon 
anomalous moment. For the one particle case, we show 
results for C v and having the same sign. Masses 
below about 160 MeV are not allowed for this scenario. 
For two particles with equal masses, the disallowed re- 
gion extends to about 210 MeV. All higher masses are 
allowed. 



V. CLOSING COMMENTS 

Exotic, in the sense of presently undiscovered, parti- 
cles that couple to muons and hadrons but not electrons 
could be responsible for an extra energy shift in muonic 
hydrogen, and thereby reconcile the muonic and elec- 
tronic proton radius measurements. However, for exotic 
explanations to work, there are requirements to be met, 
a number of which have been discussed already. 

So far we have not confronted all the additional con- 
straints placed by neutron scattering and meson decays 
mentioned in Q . To note some of these briefly, there are 
decay constraints following from searches for unknown 




Af* (MeV) 

FIG. 6: Mass limits on polar and axial vector particles due 
to constraints placed by K — > fiX searches. The solid curve 
is the result for a single particle with both polar and axial 
vector couplings, accounting for the experimental efficiency, 
with couplings obtained through satisfying the Lamb shift and 
magnetic moment criteria. The dashed curve is the result for 
separate polar and axial vector particles with equal masses. 
The experimental limit is the horizontal line, and the shaded 
region is allowed. 



particles in T or J/ip decays. These are potentially seri- 
ous. However, they do not apply if there is no coupling 
to heavier quarks. Hence a new force or new particle 
that spoofs a smaller proton radius in muonic hydrogen 
should couple only or almost only to second generation 
leptons and first generation hadrons. In particular, there 
can be coupling to muons but not electrons and to first 
generation quarks but not b or c quarks. 

Further noted in Q, neutron scattering constraints 
only limit models with very light new particle masses 
(under 5 MeV) , other muonic atom energy splittings give 
bounds already below ones discussed, the tt —> 7V decay 
where V is a massive vector only impacts an my mass 
range that is already excluded, and the limits from non- 
observance of the decay 77 — > VV is not serious for my 
below m v /2 given the results in Figs. Q] and [2] 

Fine tuning is needed to evade constraints from the 
muon magnetic moment measurements. Tuning is possi- 
ble because polar vector and axial vector exchanges give 
opposite sign contributions to {g — 2)^, and the same is 
true of scalar and pseudoscalar exchanges. The need for 
fine tuning has the additional effect that once a vector (or 
scalar) exchange contributing to the Lamb shift requires 
a further axial (or pseudoscalar) exchange in (g — 2) M , 
the axial vector (or pseudoscalar) will then potentially 
contribute to other processes, adding for example to the 
rate that needs to be evaded in a decay process. 

In models where the new physics proton and muon 
couplings have similar magnitudes, if-decay constraints 
must be avoided by having the new particle mass above 
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at least 150 MeV for the polar vector-axial vector equal 
mass case, or by having mass low, below 100 MeV or 
high, above 180 MeV, in the scalar-pseudoscalar case. 

Relatively new, much from the past year, are limits 
arising from searches for dark photons, which are massive 
vector particles that couple to charged par ticles at a re- 
duced rate. The combination of KLOE [lj|, APEX 
Mainz Al [2(|, and BABAR [2l| limit a new vector cou- 
pling to below about a few x 10 -3 of the normal photon 
coupling for the my mass range 60 to above 500 MeV, 
and below 10 -3 over part of that range. This is relevant 
to eliminating proposed theories where the magnitudes of 
the electron and proton couplings are the same or simi- 
lar. Then to give a muonic hydrogen energy sufficient to 
spoof the proton radius result would require a muon cou- 
pling so large as to violate the iGdecay bound by more 
than an order of magnitude, at least in the mass range 
60 to 310 MeV (where the phase space for the relevant 
experiment runs out). 

Our present work should be viewed proof of con- 



cept rather than a completed model. For successful ex- 
otic explanations of the proton radius problem, there are 
requirements of coupling only to targeted leptons and 
hadrons, of hne tuning the muon couplings, and of re- 
stricted mass ranges to avoid conflict with unobserved de- 
cays. While these requirements may seem difficult, there 
is still a window of possibility for new physics explana- 
tions of the proton radius problem. One would like an 
ab initio new physics theory that works for this problem, 
and the restrictions given here may be a help in finding 
one. 
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